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Abstract 

We derive spin operator matrix elements between general eigenstates of the superin- 
tegrable ZTv-symmetric chiral Potts quantum chain of finite length. Our starting point 
is the extended Onsager algebra recently proposed by R.Baxter. For each pair of spaces 
(Onsager sectors) of the irreducible representations of the Onsager algebra, we calculate 
the spin matrix elements between the eigenstates of the Hamiltonian of the quantum chain 
in factorized form, up to an overall scalar factor. This factor is known for the ground 
state Onsager sectors. For the matrix elements between the ground states of these sectors 
we perform the thermodynamic limit and obtain the formula for the order parameters. 
For the Ising quantum chain in a transverse field (A^ = 2 case) the factorized form for the 
matrix elements coincides with the corresponding expressions obtained recently by the 
Separation of Variables Method. 

Keywords: Quantum integrable spin chain, Order parameter, Onsager algebra 



1 Introduction 

The solution of the Ising model by Onsager [1] in 1944 was a major breakthrough in 
mathematical physics. The subsequent calculation of the Ising magnetization by Yang [2] 
in 1952 (after the announcement of the formula by Onsager [3] in 1949) was an achievement 
of similar importance. The main tool in Onsager's work [1] was the introduction of the 
algebra which now carries his name. Later many alternative methods for solving the Ising 
model were invented and the interest in Onsager's algebra faded away. It was not revived 
until 1985 when the superintegrable chiral Potts quantum chain (SCPC) was formulated 
[4]. For more than 40 years the Ising model had remained the only known representation 
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of Onsager's algebra. Then the SCPC provided a whole series of new Onsager algebra 
representations. In many respects, the SCPC is the TLjq generalization of the Ising 
model. Soon, in 1988, Baxter [5], using functional relations of the 2-dimensional integrable 
chiral Potts model [6, 7] calculated the partition function of the SCPC. The formula for 
the order parameter (the generalization of the famous Onsager-Yang result for the Ising 
model) was conjectured in 1989 [8] from a perturbative calculation up to k for arbitrary 
N (see also [31]). However, it took until 2005 that Baxter [10, 11] succeeded to prove the 
conjectured formula. His derivation used functional relations for the T2-model and natural 
analytical assumptions. 

Nevertheless, Baxter [12] insisted that a purely algebraic derivation of the order pa- 
rameter formula would be desirable. So, recently, he extended the Onsager algebra by 
a SCPC-spin operator [13]. In order to close the algebra one has to include many new 
operators which are produced by commutations with the generators of the Onsager al- 
gebra. R.Baxter conjectured that the two simplest relations in the extended Onsager 
algebra uniquely define (up to a common multiplier) the matrix elements of the spin op- 
erator between vectors from the ground state Onsager sectors. Baxter suggested that this 
conjecture will imply the determinant representation (7.11) in [12] which in the thermo- 
dynamic limit will give the formula for the order parameters. 

Our starting point for this research was that we compared the factorized formula 
(78) in [14] for the Ising quantum chain spin operator matrix elements obtained by the 
Separation of Variables Method to Baxter's conjectured equation (3.45) of [13]. In the 
case of the Ising chain, after a simple transformation given in Section 4.1 of the present 
paper, (3.45) of [13] should imply (78) of [14]. Since Baxter's conjecture was formulated 
for arbitrary parameters, it was natural to suggest that the factorized formula (78) in [14] 
for the ground Onsager sectors is valid also for the ground Onsager sectors in SCPC if one 
uses appropriate parameters labelling these sectors. Numerical calculations for the finite 
length chain together with calculations for the thermodynamic limit (see Appendix C) 
strongly supported our conjecture. This encouraged us to prove the Baxter's conjecture 
(3.45) from [13]. 

In the present paper we generalize this conjecture to arbitrary Onsager sectors and 
prove it. We found that the matrix elements of spin operators between the eigenvectors 
of the SCPC Hamiltonian for arbitrary Onsager sectors can be presented in a factorized 
form. The problem of finding the matrix elements of the spin operator is reduced to the 
problem of finding non-zero matrix elements of the spin operator between some particular 
vectors from the corresponding two irreducible representations of the Onsager algebra. 
In the case of ground state Onsager sectors this gives an exact factorized formula for 
the matrix elements of spin operators between such type eigenvectors. Finally, taking 
the thermodynamic limit of ground state matrix elements we derive the Albertini et al. 
formula [8, 31] for the order parameters. We would like to note that factorized formulas 
for the spin matrix elements exist also for the 2D Ising model [15] and for the quantum 
Xy-chain [16]. 

Recently, state vectors for the superintegrable chiral Potts model were investigated 
exploiting the s/(2)-loop algebra symmetry of the T2-model [28, 29], which is closely related 
to the chiral Potts model [17, 18, 19, 20]. Since our approach avoids direct use of the 
relation to the T2-model, we shall not go into details of this exciting work. 

This paper is organized as follows: In Section 2 we define the SCPC by its Hamiltonian 
and use the fact that its two parts generate the Onsager algebra and imply the Ising form 
of the eigenvalues. Baxter's extension and the basic two resulting equations for the spin 
operator are stated. In Section 3 we develop the consequences of Baxter's two equations 
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for spin matrix elements: including higher Onsager sectors and in several steps fixing the 
explicit solution. Specializing to the ground state sector, we show that our result coincides 
with Baxter's conjecture. The Section ends with the discussion of some non-trivial selec- 
tion rules for the spin matrix elements. In Section 4 we sandwich the expression for the 
spin operator between eigenvectors of the Hamiltonian. The sum over the intermediate 
labels is performed explicitly, resulting in fully factorized expressions. Finally, performing 
the thermodynamic limit in Section 5, the formula for the order parameters is obtained. 
Section 6 summarizes our results. Technical details about performing the factorization 
and the thermodynamic limit are given in Appendices, together with a numerical example. 



2 The superintegrable chiral Potts quantum chain 

2.1 The Hamiltonian of the superintegrable chiral Potts 
quantum chain 

The superintegrable ZAr-symmetric chiral Potts quantum chain of length L is defined by 
the Hamiltonian [4] 

n = nQ + v!'Hi, (1) 

where k' is a temperature-like real parameter, oj = e^'^^l^ ^ A'^ G Z, N >2 , and 

L N~l yn y—n L N—1 v-n 

j=l n=l j=l n=l 

The space of states is the L-fold tensor product of spaces at each site j. The operators 
Zj, Xj G End(C^)'^'^ act non-trivially only at site j: 

Zj = 1 (S) ...^ 1 (g)... (g) 1, 

j-th 

= 1 . . . (g) 1 (g^A^^l 0... (g) 1. 

j-th 

Z and X are the Zjv-generalizations of the Z2-Pauli operators az and ax ■ They satisfy the 
Weyl relation Z X = uj X Z and Z^ = X^ = 1. We always take periodic boundary 
conditions Zl+i = Zi. A convenient representation is to label the state vectors by spin 
variables aj G Z^v, j = 1 ■ . . L, writing the state vectors as |cri, CJ2, . . . , CTl) and 

(^j)<7j,a; = ^cj.^a'.^"' , = K,,a'^+1 ■ (3) 

Both Ho and Tli commute with the spin rotation operator R = Xi X2 ■ ■ ■ Xl- Since 
= 1, the space of states of SCPC decomposes into sectors of fixed of ZAr-charge Q 
where Q = 0, 1,...,A^ — 1, corresponding to the eigenvalues uj^ of the R. 

It turns out that < k' < 0.902... is the ferromagnetic regime [33], where in the 
thermodynamic limit the system acquires a non-zero value of the order parameters. One 
aim of this paper is to derive its value as a function of k'. 

The Hamiltonian (1) can be derived as a logarithmic derivative with respect to a ra- 
pidity parameter from the homogeneous 2D superintegrable chiral Potts model. However, 
we shall avoid using properties of the two-dimensional model. 
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2.2 The Onsager algebra and the superintegrable chiral Potts 
quantum chain 

In the beginning of the 1980'th chiral Potts quantum chains received much attention since 
for certain parameter values they exhibit an incommensurate phase. Since no integrable 
chain of this class was known, Howes, Kadanoff and den Nijs [21] studied perturbative 
expansions. They found strong evidence that the Z3 version of (1) has a terminating 
low-temperature expansion for the mass gap. Subsequent numerical studies in [4] showed: 
if in the ZAr-symmetric quantum chain the coefficients of the terms involving Zj Zj^i and 
Xj are chosen as in (2), many gaps become linear. This feature is known from the Ising 
model and motivated further studies of the Hamiltonian (1), (2). Now is was known that 
Hamiltonians of the form 7^ ~ Aq + k' Ai are integrable if Aq and Ai satisfy the 
Dolan-Grady-relations [22] 

[Aj, [A„ [A„ ]]] = 16 [Aj,A^^j ] , j = 0, 1. (4) 

Indeed, in [4] it was shown that, putting 

Ao = -2ni/N , Ai = 2no/N . 

the relations (4) are satisfied. A pair of operators Aq, Ai satisfying (4) generates a basis of 
Onsager's infinite dimensional Lie-algebra with the operators Am, m G Z, Gn, n G N : 

[^rrt) ^n] — ^Gm—n i 
mi -^n\ — ^Am+n "^Afi.—m , \Guii Gn\ — . (5) 

The study of the representation theory of the Onsager algebra has been started by 
B. Davies [23]. He showed that certain finite-dimensional representations of the Onsager 
algebra are related to representations of a direct sum of several copies of the Lie algebra 
sl{2), see also [24, 25, 26]. In the case of the SCPC only two-dimensional representations 
of sl{2) appear. Unlike the N = 2 Ising case, for > 3 the Onsager algebra alone is not 
powerful enough to fully determine the eigenvalues. In each Onsager sector of the SCPC 
the 2*"^ eigenvalues of Ti have the Ising-like form 

rriE 

E = A + V!B - N^±e{ej) , (6) 
i=i 

with 

e(^) = Vl - 2k' cos + k'2 . (7) 

The constants A and B and the niE values 9j label irreducible representations of the 
Onsager algebra. E.g. for = 3 values of niE can be found in [27]. However, detailed 
formulas for niE will not be used in the following. 

Using functional relations which follow from the general integrable 2-dimensional chiral 
Potts model [29], or, recently, exploiting the loop symmetry of the related T2-model [17], 
the 9j can be calculated from the zeros of the Baxter-Albertini-McCoy-Perk (BAMP) 
polynomials [5, 30, 9]. However, here we shall follow the spirit of Baxter's recent work and 
determine the spin-operator matrix elements from the Onsager algebra and an extension 
of this algebra alone. So, the 9j will be considered as free parameters. Of course, specific 
values of the 9j calculated from the zeros of the BAMP-polynomials may be inserted in 
the final formulas. 
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2.3 Extension of the Onsager algebra by a spin operator 

Baxter proposed [13] to extend the Onsager algebra of the SCPC by the local spin operator 
S = where r may be chosen to take the values r = 1, . . . — 1. He introduced two 
operations on operators X in the space of states 

MX)^^^. and /,,X,.£^L±^ (8) 

and has shown that 

/o(5) = 0, and fiifiiS)) = fi{S) . (9) 

Repeated commutation of S with Tio and 7ii produces many new elements. However, 
Baxter made the conjecture that in the case when the finite-dimensional representations 
of Onsager algebra are related to two-dimensional representations of sl{2) (in particular, 
in the case of the SCPC) already the mentioned two relations for S uniquely (up to a 
multiplier because the relations are homogenous) define the matrix elements of S. In the 
following section we prove this conjecture. 

3 Conjecture on the reduced matrix Spq: a gen- 
eralization, reformulation and a proof 

3.1 Generalization 

We consider the restriction of Hamiltonian Ti on the subspaces which are spaces of ir- 
reducible representations of the Onsager algebra. The Hamiltonian 7i on these spaces 
has an Ising-like spectrum. We call these spaces Onsager sectors. Let , be some 

p Q 

Onsager sectors with Zjy-charges P, Q and dimensions 2'"b, 2'"b, respectively. We will 
use the short notations m = rn^, n = m^. 

We denote r = Q- Pioz:P<Q, and r = Q- P + NioTP>Q. For the sector 
we have 

HP = H^ + k'H[, H^ = AP-NF^, H^ = bP + NF[, 
where and are some constants depending on this Onsager sector and 

" ^ f -1 ) , ' \ sin^,- -COS. J 



We write the same formulas for the sector with charge Q and with parameters Q'^^ j 
1,2, ... ,n: 



n 



^0 Z^l -1 / .' 1 sine' -cos^' , . 



Using the fact that Hq and Hi are block-diagonal with respect to irreducible representa- 
tions of the Onsager algebra and taking the formulas (8) between arbitrary states from 
and V*^, we obtain 

/o(Xpq) = -hp^XpQ - XpqF^ + aXpQ) , (10) 
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MXpq) = ^{FfXpQ - XpqF? + {(3 + 1)Xpq) , (11) 

where Xpg is rectangular submatrix in X and 

a = {A'^ -A^)/N, f3 = {B^ -B^ + 2r-N)/N . (12) 

Therefore it follows from the relations (9) that the rectangular submatrix Spq of the size 
2™ X 2" of spin operator S = Z\ satisfies the system of equations 

F^C - CF^ + aC = , (13) 

Ff Ff C - 2F[CF^ + CF^F^ + 2(3{F[C - CF^) + {(3^ - 1)^ = (14) 

with respect to the unknown matrix C. The F^^, Fq , Ff, F^ are given above and a, j3 
are supposed to be arbitrary constants. Formulas (10), (11), (13), (14) generalize Baxter's 
relations from [13] which were written only for ground state Onsager sectors. 

The system of equations (13), (14) is homogeneous and any solution of this system 
multiplied by a constant is also a solution. We will find a solution Cpq of the system 
(13), (14) and show that space of solutions is no more than one-dimensional. Therefore 
exact^ matrix elements Spq of the spin operator S = Zl are related to this solution by 
relation Spq = Mpq Cpq with an unknown scalar multiplier Afpq, which should be fixed 
from other arguments. The other relations of the extended Onsager algebra do not spoil 
the equality Spg = MpqCpq, but they may impose that Mpq = 0. 

3.2 Reformulation 

Multiplying the relations (13), (14) from the right by 

n 

G = J{<T-'y-, (15) 
i=i 

using Fq = —G~^FqG, F^ = —G~^F^G and denoting C = CG, we get similar relations 
but with changed signs at Fq and F^: 

F^C + CF^ + aC = , (16) 

F[f[C + 2F[CF^ + CF^F^ + 2(3{F[C + CF^) + {(3^ - l)C = . (17) 

Let us combine the matrix elements C^,/? into a vector K S ® with components 
labeled by multi- indices. The vector K has dimension 2^, where fi = m + n. Let 
for 'k = 0,l 

Ffc = Ff ® 1 + 1 (Ff )* 

be an operator in and {F^Y is the transposition of F^ . Then the relations (16) 
and (17) become 

{Fo + a)K = 0, {Fi+P-l){Fi+P + l)K = 0. 

This procedure puts the parameters {9j} and {9j} on an equal footing: we combine these 
two sets into the set . . . , t?^} = {^i, . . . , 6m, O'n ■ ■ ■ , O'n}- 

^In the following the term "exact matrix elements" will be used in the sense that the constants Afpg have 
been determined too. 
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Since the matrices Fq and Fi are both traceless with determinant —1 we can find a 
similarity transformation Fi = Fq U to get 



where 



{Fo + P-l){Fo + p + l)UK = 0; 



cosi9j + 1 cos??j — 1 \ TT — TT f ^ ^ 



(18) 

It gives finally two relations for the vector K = UdK: 

{Fo + a)K = 0, {Fo + (3-l){Fo + /3 + l)UK = 0, (19) 

where we use the fact that -Fq commutes with Ud- In the next two subsections we will 
solve the system of equations (19) with respect to K for arbitrary parameters a, (3. 

3.3 Number of solutions of (19) 

The relations (19) are a system of linear homogeneous equations with respect to the 
unknown components of K. In this subsection we estimate the number of solutions of this 
system, depending on a and (3. We use the notation F = Fq. 
Theorem 1. The dimension of the space of solutions K of the system 

{F + a)K = 0, (20) 

{F + p-l){F + P + l)UK = 0, (21) 

where 

j=l ^ ^ J j=l^ ^3 

at generic values of the parameters c^ = cos'&k, k = l,...,/z, is no more than one- 
dimensional. It is one- dimensional if and only if 

• ji is even, a = 0, (3 = ±1, (23) 

• 12 is odd, a = ±l, /? = 0. (24) 



Proof. The spectrum of F is fi, fi — 2, . . . , —fi. From the first relation (20) we have a 
necessary restriction on a in order to have K ^ 0: a should be an integer of the same 
parity as n, and \a\ < fj,. If a satisfies this condition, the space of solutions of (20) has 

dimension ( , ^ n ,„ I • The second relation (21) requires the existence of such K from 
V (^ + a)/2 ; y J ^ 

this subspace which are annihilated by the operator {F + (3 — 1){F + (3 + 1)U. Thus, to 

have K ^ 0, ^ + /3 should be an odd integer and < /-f + 1. Otherwise, since U is 

invertible, we would have that {F -\- (3 — 1){F -\- (3 -\- 1)U is also invertible and therefore 

K = 0. 

We will prove the theorem by induction over fi. Let F = Fj-^) be the initial matrix. In 
a basis in which the ^-th tensor component is separated, we have 



1 \ / F(^_i) + Q + 1 



(,)^a = .(,_i)^ -1 F 



+ a 
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It is a diagonal matrix and the zeros on its diagonal define the space of solutions of (20). 
Let P(^) (a) be the projector to this subspace. Its rank is 

In the basis with separated fi-th tensor component this projector is 

We need to calculate the rank 

f{fi,a,(3) = rk +(3- +(3 + l)C/(^)P(^)(a)) , 

since the difference 

fifi, a, (5) = rk (a) - /(/x, a, /3) > 

defines the dimension of the space of solutions of the system of equations (20) and (21) 
Extracting explicitly the /i-th tensor component we have 



f{fi,a,P) = rk ^ ^ ^ ^ 



where 

A = {c^ + l)(F(^-i) +P + 2)(F(^_i) + /3)C/(^_i)P(^_i)(a + 1) , 

B = {c^- l)(F(^_i) +P + 2)(F(^_i) + /3)C/(^_i)P(^_i)(a - 1) , 

C = (F(^_i) + /3)(F(^_i) +P- 2)C/(^_i)7'(^„i)(a + 1) , 

D = (F(^_i) + /?)(F(^_i) +P- 2)[/(^„i)P(^_i) (a - 1) . 

By f{fj,,a,(3) we denote the rank which is calculated at generic values of the parameters 
Cfc, A; = 1, . . . , /i. At some specific values of these parameters, the corresponding rank 
may be smaller than at generic values. At = 1 we have B = and using the fact that 
the rank of a block-triangular matrix is greater or equal to the sum of ranks of diagonal 
blocks we obtain 

/(//, a, /?) > rk A + rk D = /(;u - 1, a + 1, ;9 + 1) + /(;u - 1, a - 1, /3 - 1) . (25) 
Similarly, at = — 1 we have A = and 

fin, a, /3) > rk S + rk C = /(/i - 1, a - 1, /3 + 1) + /(/i - 1, a + 1, /3 - 1) . (26) 

Using identities between binomial coefficients, the inequalities (25) and (26) become 

f{li,a,(3) < f{i^-l,a + l,(3 + l) + f{i^-l,a-l,(3-l), (27) 

/(/^,a,/3) </(/U-l,a-l,/3 + l) + /(;U-l,a + l,/?-l). (28) 
We want to show by induction in ji that 

/(^, a, /3) < .5,2+^2,1. (29) 

After proving this inequality for fi — 1 and all a, f3, a straightforward analysis shows that 
the inequalities (27), (28) imply (29). 
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Assuming that the binomial coefficients ^ ^ J are defined to be if Z < or / > /c, 

we can start the induction from the case /i = 1 and arbitrary q, (3. In this case, as it has 
been explained in the beginning of the proof, in order to have X 7^ we need to choose 
a = ±1 and (3 S {—2, 0, 2}. Let for definiteness a = 1. Then 

/(l,a,/3) =rk ((F(i)+/3-l)(F(i)+/3 + l)C/(i)P(i)(a)) = 
_ / (/3 + 2)/3 \rj( ^ ^\ 

\ - 2) ; ^ 1^ 1 ; ■ 

If /5 = then /(l,a, /3) = 0, otherwise /(l,a,/3) = 1. Thus only if /3 = we have 
/(l,a,/J) = l. 

This proves that the space of solutions K is no more than one-dimensional if a and (3 
satisfy the relations (23) or (24). Otherwise there is only the trivial solution K = 0. In 
the next subsection we give the explicit solution for K \i a and /3 satisfy (23) or (24). It 
proves that in these cases the space of solutions is one-dimensional. □ 

What is the condition on the parameters c^, k = 1, . . . which leads the theorem 
to fail? The analysis for < 3 gives that this condition is: either two of Cfc coincide and 
equal to ±1 or three of Cfc are pairwise coinciding. We conjecture that for general /i, the 
corresponding condition is too restrictive and that the sets arising from SCPC do not 
satisfy this condition. 



3.4 The solution for the matrix elements 

We proved that if a non-trivial solution K exists, then it is unique up to a multiplier. 
Such a solution exists if and only if is even, a = 0, /3 = ±1, or if /i is odd, a = ±1, 
/? = 0. 

We now present the solution K in the basis of eigenvectors of F. Each basis vector 
is a tensor product of the two-dimensional vectors associated with summands of F in 
(22). We label each basis vector by the subset I of the set {1,2,...,//}, such that the 
elements of I label the tensor components where the two-dimensional basis vectors with 
the eigenvalue +1 are located. The complementary subset / contains the labels of the 
tensor components where the two-dimensional basis vectors with the eigenvalue —1 are 
placed. By |/| we denote the size of the subset /. 

We use a short notation: the product of an expression containing the variable c indexed 
by a subset / means the product of this expression over the elements of the subset with 
variable c indexed by the elements of /. 

Theorem 2. Let the parameters a and (3 satisfy (23) or (24)- The vector K with the 
components 

^ ^ V|,(M-")/2 /onN 

' \{{ci-ci)\[{ci + lY\{{ci-lY ^ ^ 

is the unique (up to a scalar multiplier) solution of the system (20) and (21). Here 

(J = (/? - a + l)/2, T = {(3 + a+l)/2. (31) 



Remark. Due to the special values of a and /?, each of a and r can be either or 1. 
Proof. The numerator '^|/|,(^_o)/2 ensures that K satisfies (20). 
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Let us prove that K satisfies (21). In fact we need to prove that the non-zero compo- 
nents of M = UK are in the eigenspaces of F with eigenvalues —(3 — 1 and — /? + 1. The 
components of the vector M = UK are 

M \^rr k- ^ ^(c'nj + 1) rK^nj - l)i'|/U;.-a)/2 , , 

M. = ^ u,jK, = ^ n(<:.-c,-)nfo + i)" n(cr-i)' ■ 

where we used Ujj = Yl{cjnj + 1) Yl{cjf^j — l) which fohows from (22). The components 
of M are 

P 

^ Uicj-cj) uicj+ir Uicj-iy ^^^^ 

with an unknown polynomial P in the numerator. Let us explain the origin of the different 
factors in the denominator of (33). First, there are no factors like Cj + 1, j G J in the 
denominator of (33), because all such factors with j & I in the denominator of (32) 
are cancelled. Similarly, there are no factors like cj — 1, j € J. Second, there are no 
factors like cj — Cj^, € J, in the denominator of (33). Potentially such factors can 
arise from summands in (32) for which j € / and j' € /. But to each such summand there 
corresponds a summand with / replaced by / = {j'} U I\{j}. A straightforward analysis 
shows that the sum of these two summands does not contain a pole at cj = Cji. Similarly 
there are no factors like Cj — Cj', € J in the denominator of (33). 

Let us estimate the degree of P with respect to the variables Cj, j € J. It is 

max(l - a - |/|, 1 - r - |/|) + |J| > degP > 0. 

The two arguments of the maximum correspond to the fact that each element of J is 
either in / or in /. The term \J\ corresponds to the denominator of (33). Similarly, the 
estimate of the degree of P with respect to the variables Cj, j G J, gives 

max(-cr- \I\,-T - \I\) + \J\ + a + T > degP > 0. 

In order to get M = UK ^ both inequalities have to be satisfied. Using |/| = {fi — a)/2, 
\A = (^ + Q^)/2, these give 

\J\>{^i + P-l)/2, |j| > (^-/3_i)/2. 

Since |J| + |J| =// there are only two possibilities: |J| = (/x — /3 — l)/2, |J| = (;U + /3 + l)/2 
or \J\ = {fi — (3 + l)/2, \ J\ = {fj, + f3 — l)/2. The corresponding eigenvalue | J| — | J| of F 
is — /? — 1 or —(3+ 1, respectively. It means that K satisfies (21). The theorem is proved. 
□ 

In order to obtain matrix elements of 5 = Z[ between the states from Onsager sectors 
we need to rewrite the solution K given by (30) in the original notations. 

We will label each basis vector in by the subset V of the set {1, 2, . . . , m} cor- 
responding to the tensor components where the two-dimensional basis vectors with the 
eigenvalue +1 of are located. The complementary subset W contains the labels of the 
tensor components for which the two-dimensional basis vectors with the eigenvalue — 1 
are placed. In the same way we define the subsets V and W of the set {1, 2, . . . , n} to 
label the basis vectors in the space . 

We also use the set of indices {1,2,..., m+n}. We will identify the indices {1,2,..., m} 
with indices of tensor components of . The rest indices {m + 1, m + 2, . . . , m + n} of 
the set {1, 2, . . . , m + n} are identified with indices {1, 2, . . . , n} of tensor components of 
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by subtraction m. This identification defines a one to one correspondence between 
the subsets in {1, 2, . . . , m + n} with pairs of subsets from the sets {1, 2, . . . , m} and 
{1,2,... 

Using the relation C = CG with G defined in (15) and the relation K = U^K with Ud 
defined in (18) we have for the matrix element of G in the basis labeled by V and V' 

Cvy = (-1)'^ ^Cvy = (-1)'^ '^(y,iy) 

= (-l)l^li<:/ JJ tan 0^/2 JJ tan 0^2, (34) 

where / = {V, W') is the subset in {1, 2, . . . , m + n} defined by the subsets V and W' 
according to the identification given above. 

Note once more, that Theorem 2 defines the matrix elements {GpQ)vy by (34), (30) 
as solutions of the system of homogenous equations (13) and (14). Exact matrix elements 
{SpQ)vy of the spin operator are given through the relation 

{SpQ)vy = ■N'pQ{GpQ)vy (35) 
with an unknown constant Mpq which depends on both Onsager sectors. 



3.5 Matrix elements between ground state Onsager sectors 

For the ground state Onsager sectors we have 

= L{1- N)+Nm, A^ = L{1- N) + Nn, (36) 

= 2P + L{l-N) + Nm, B'^ = 2Q + L{1 - N) + Nn . (37) 

It was argued in [13] that in the ground state Onsager sectors one has Syyi = 1 for 
= 0, W' = 0. In particular this relation defines exact values of the constants Mpg for 

all pairs of ground state Onsager sectors. 

Let us consider the case P < Q, m = n. Using (36), (37), (12) and (31) we obtain 

II = 2m, a = 0, f3 = —1, a = t = 0. The components (30) of K are 

Ki = TTp^ , (38) 

n(c/ - cf) 

where c/ is the set {cos 6*4, cos 6*^}, i G U, j € W and cj is the set {cos cos 0^}, i G W, 
j € v. Now using (34) and (38) we get Gyy- From the relations Syy = GyyMpQ and 
Syy, = 1 for = 0, VF' = we derive 

■^"^ " nr=itan^;./2 • 
Thus our solution for Gyyi implies 

Svy = ^\w\,\W'\ tan6'j/2 cot6'-/2 



T{i(iV,j&wi^i ~ ^j) T{i(^W'j(^V'i'^i ~ 
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This formula for Syy coincides with the formula (3.45) conjectured by R. Baxter in 
[13] for the case P < Q, m = n. It means that Theorem 1 and Theorem 2 prove this 
conjecture. 

Let us consider the case P < Q, m = n + 1. We have ji = 2m — 1, a = —1, /3 = 0, 
cr = 1, T = 0. In this case the components of K are 

Kr = . (40) 

Now using (34) and (40) we get Cyy = SvyMpq, where 

Svy = \w\,\W'\ W sin6'i (sin6'-)"^ 



,1 \ ' 



This formula for Syyi is again Baxter's formula in the case P < m = n + 1. 

We have also two other cases corresponding to P > Q which can be considered analo- 
gously. In these cases our calculations also prove the Baxter's formula for Syyi. The first 
case is P > Q, m = n and we have fi = 2m, a = 0, P = 1, a = t = 1, 

m n m n 

NPQ = {-irj{{l + c.) n^in^;- \{J{{c^-c'^)■ (42) 

i=l j=l i=l j=l 

The second case \s P > m = n — 1 and we have /i = 2?! — 1, a = 1, /3 = 0, o" = 0, r = 1, 

n m n 

AfpQ = i-irUsine'^ nn(^^-^^)- (43) 
j=l i=l j=l 

3.6 Application: selection rules for spin matrix elements 

We proved that the non-trivial solution exists if and only if fj, is even, a = 0, (3 = ±1 , or 
if /i is odd, a = ±1, (3 = 0. This means that only in these cases we may have non-zero 
matrix elements Spq = MpqCpq for spin operators. These selection rules are more fine 
than the selection rule by Z^v-charge. They follow from (12). For r = Q — P for P < Q, 
and r = Q — P + N for P > Q we have two cases where Cpq ^ 0: 

• ^ = m + n is even and - = 0, - B'^ = N ± N - 2r. 

• fi = m + n is odd and A*^ - A^ = ±N, B^ - B^ = N - 2r. 

To be precise, these selection rules are not complete because the matrix elements for 
spin operators between vectors from two Onsager sectors were found up to a common 
multiplier Npq. For a pair of specific Onsager sectors this multiplier can be zero even if 
CpQ 7^ 0. Additional information is needed whether Mpq is zero or not. 

We have verified these selection rules numerically for N = 3, L = 3. In this case they 
are complete, that is all Mpq 7^ for the Onsager sectors defined by these selection rules. 
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4 Factorized formula for matrix element of spin 
operators for a finite chain 

4.1 Matrix elements of spin operators between the eigen- 
vectors of the Hamiltonian 

As has been explained in Section (3.1), the Hamiltonian of the model restricted to the 
sector with charge P is 

H^ = A^ + k'B^ -Ny( ^- ^- r^'T^\ 

^\ -k'smOj -l + k'cos6'j 



The eigenvalues of are 



E = AP + k'B^ -N^±e{ej) 



where 



£{9) = Vl - 2k' cos + k'2 , e(0) = l-k', e(7r) = 1 + k' . 



The —e{6j) (resp. +e{9j)) in the sum corresponds to the presence (resp. absence) of the 
excitation with rapidity 9j because the change e{6j) to —s{Oj) increases the energy. The 
eigenvector of with the lowest energy (that is without excitations) is 

ai(0i) ai(02) ® V 



ai{e) 



k' sin 61 



k'2 sm'^e + {l-W cos9-e{9)f 

1 - k'cos6'-e(6i) 
k'2sin2^ + (1 - k' cos^-e(0))2 

al{e) + alie) = i. 



It is useful to express the components ak{Oj), k = 0,1, of the eigenvector through the 
energies e{6): 



ai{e) 



/ (g(g)-g(0))(e(7r)-e(g)) 
Ae{9) 



The formulas for the eigenvectors with some 9j excited are obtained from the formula 
for the non-excited eigenvector by replacement £{0j) — > —e{9j) for all such Oj. Therefore 
it is sufficient to consider matrix elements of spin operator only between non-excited 
eigenvectors from each Onsager sector and to make the replacement £{0j) — > —£{9j) for 
all excited 9j in the final formula. We use the same formulas for the sector with charge 
Q and with parameters 9j, j = 1,2, . . . , n. 
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In each Onsager sector we call the eigenvector of H with the lowest energy (that is 
without excitations) the ground state and denote it \g.s.)p. The matrix element between 
the ground states of two Onsager sectors with charges P and Q, r = (Q — P) mod A^, is 

p{g.s.\Z{\g.s.)Q = Y.ll J] J] «o(^D H '^MiSpQhy • 

V,V' i&V ieW i£V' i£W' 

The exact matrix elements {SpQ)vyi are related to the matrix elements {CpQ)vyi defined 
by (34), (30) by the relation (35). Then 

p{9-S-\Zl\g.s.)Q 

^ M,,f[Me.)Y[am E n (-^-4) n (^-4) 

i=i 1=1 vy'iGW \ uv / -gy, \ i\ ij / 

= (_i)" _/^^^(2k')''^^""^/^+^^^+"^/^+^''' f[ /^i£i±£oll£zL±£il^ x 

i=l ^ ^ 



I i€l 



where we defined = e{9i) for 1 < i < m, = —e{9^_^) for m < i < m + n and 
eo = e(0), = £(^)- Also we used 



cos 0- cos 0' = {e\9')-e\9))/{2k'), tan 



9_ £^{9)-e^0) 
2~\1 e^{TT)-e\9) 



9 ai{9) _e{9)-e{^) 9 aQ{9) _ -e{9) - e{{)) 

^^"^ 2 ■ ~ e(7r) + e{9) ' *''''2'^r(^" ^(vr) - e(0) ' 

In order to perform the summation over I in (44) we have to discuss the cases of even 
and odd separately. 

4.2 Summation in the case of even fi 

In this case we have a = 0, P = ±1. Let us consider the case /3 = —1. We have a = t = 0. 
We use the following summation formula: 

^UM^ ^ 

(_1)A'(M-2)/8(^^ _ 
YlciZc + u){Zc + V) UcKsi^s + Zc) ' 

which is proved in Appendix A. This formula at Zi = Si, i = 1, jj,, and u = —Sq, f = 
gives 



/ Uieii^i - ^o) Uieii^^ + ^i) UieidGii^] - 4) 

(_l)MM-2)/8(g^^g^)/./2 

nf=ife - ^o)(e^ + Si) Ui<ji^i + ^j) ' 
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which we use to make the summation in (44): 
p{9-S-\Zl\g.s.)Q = 



m 



X 



(_l){''n+n){m+n-2)/8+n(n+l)/2^2k')('"+")^/^ 

Now let us consider the case of matrix elements between ground states of the ground 
state Onsager sectors for which P < Q and m = n = ijl/2. In this case we can use the 
expression for Mpq from (39) rewritten in terms of e{9) 



= (-i)-(2ko— n \l e^(e'^)-eko) ■ n n(^(^')' - ^(^^•)') 



to get exact matrix elements of spin operator between ground states 



In the case of general m, n and a = 0, /3=1, (t = t = 1 the analogue of (45) is 

m 

p{9-S-\Zl\g.s.)Q =7Vpq11 



(46) 



X 



\ (e(7r) - e{ei))^2e{ei){e{ei) + e(0))(£(7r) + e(0,)) 

TT ^ 

y (.(vr) + e{e[))^2e{e[){e{e[) - s(0))(e(7r) - e{e[)) 

( _-\\{m+n){m+n-2)/ii+n{n+l)/2 ( iy[,l-\m+n+{m+n)'^ / i 

(47) 



nr<,(e(^.) +^(%))nr<,(^(^D +^(^;-))n.=in-=i(e(^.) -^(^P) ■ 

This gives the matrix elements between ground states of the ground state Onsager sectors 
for which P > Q and m = n = /x/2if one uses the corresponding A/pg from (42). The 
final formula is analogous to (46) with additional factor (—1)"^ and with the replacement 
of the sets {£{ei)} ^ Me'j)}. 

4.3 Summation in the case of odd fi 

In the case of odd /x we have a = ±1, /3 = 0. Let us consider the case a = —1. Then 
o- = 1, T = 0. 

In this case we will use another summation formula from Appendix A: 



E 



1^I,(m+i)/2 



(_1)(m+1)(m-1)/8(^ _^ y)(M-l)/2 
Wci^c + u){Zc + V) nc<s(^s + Zc) ' 
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It allows to perform the summation in (44) and get 



P^g.s.\^, \g.s.)Q ^^^^^ ^ ^^^^^^/^ ^^^^^ _ ^^^^^^ y 2e{em^) + e{6,)) 



n 

X 



n 



(_l)(M'-l)/8+n(n+l)/2^2k')(''+^)'/'' 



(6(^0 - em _ 

D) 

(48) 



Now let us consider the case of matrix elements between ground states of the ground 
state Onsager sectors for which P < Q and m = n + l,m = (fi + l)/2, n = (fj, — l)/2. In 
this case we can use the exact expression for Afpq from (41) rewritten in terms of £{9): 



(-i)"»(2k')-'"» n n ^«)') n n V^Ji;^ 



i=l j=l i=l j=l V ^ J 

to get the exact matrix elements between ground states: 

p{g.s.\Zl\g.s.)Q = 



fr / (£(g.)+£(0))(e(7r)+e(g.)) -A- 1 

t\ V 2e(^,) ii Y 2e(0O(^(^D + ^(0))(e(^) + e{9'^) 

In the case of general m, n and a = l, /3 = 0, (t = 0, r = l the analogue of (48) is 

^fpQ ^ 1 



X 



(e(0) +e(7r))V2 11 ^2e(0,)(e(^.) + ^(O))(£(vr)+e(0.)) 

n 

n . 

y V2e(eO(^(^0-^(o))(e(^)-^(eD) 

■^(m+n+l){m+n+3)/8+n(n+l)/2('9^/>l(m+n+l)V4 

(50) 



It gives the matrix elements between ground states of the ground state Onsager sectors 
for which P > Q and m = n — 1 if one uses the corresponding Mpq from (43) . The final 
formula is analogous to (49) with an additional factor (—1)" and with the replacement of 
the sets {e{ei)} ^ {e{9'j)}. 
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5 Order parameters for the superintegrable chiral 
Potts quantum chain 

5.1 Thermodynamic limit of ground state matrix elements 

For the superintegrable chiral Potts quantum chain with Hamiltonian (1) we have — 1 
order parameters Sr, r = 1,2, N — 1, which take non-zero values in the ferromagnetic 
phase < k' < 1 and are defined by the matrix elements between ground states of ground 
state Onsager sectors (labeled by charge P) in the thermodynamic limit: 

Sr= lim p{g.s.\Zl\g.s.)p+r ■ (51) 

L— >oo 

At < k' < 1 there is A^-fold asymptotic (L ^ oo) degeneration of energies of these 
ground states. (In fact, for 0.901292 . . . < k' < 1 these translationally invariant states 
are not the true ground states, i.e. states with lowest energies [33, 30]. Due to level 
crossing, the lowest energy state becomes non-translationally invariant). In the paper [8], 
from perturbative calculations the formula for the order parameters 

Sr = {l- k'2)K^-r)/(2JV2) (52) 

was conjectured. Recently this formula was proved by Baxter [10, 11] using functional 
relations and analytical properties. Below we derive this formula starting from the exact 
results for the matrix elements of spin operator between ground states in the chain of 
finite length L. 

From (46) the square of ground states matrix elements in the case of m = n is 



p{9-s-\Zl\g.s.)p_^ 



nr,;=i(g(^o+^(^i))' ^^^^ 



where the 6i and 9'^ refer to the Onsager sectors with charges P and P + r, respectively. 

In the case m = n + 1 from (49) the matrix element of the spin operator Z[ between 
the ground states is 

' ''"'^^ (eiO) + e(7r)) 1^=1^) + e(0))(e(eO + e(^)) 

From (51) and (52) the formula for the order parameters which we have to obtain in 
the thermodynamic limit is 

hm p{g.s.\Zl\g.s.fp^, = (1 - k'2)^(^-^)/^' . (55) 

L—KyD 

First we consider the case of matrix elements of type (53) when m = n. Whereas 
up to now the 6i and 9'^ were arbitrary, now we take them to be determined by the 
BAMP-polynomials of [5, 30]. We use the following relation for the thermodynamic limit: 

^ log(A + e{9.)) - ^ log(A + g(gD) = ^ log ^ 1 ^ > (56) 

1=1 i=l / 
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where A is arbitrary parameter. The proof of this formula is given in Appendix B. The 
thermodynamic hmit is taken in such a way to keep the relation m = n and values of P 
and r. This limit can be realized by adding multiples of N to the chain length L. 
Using (56) at A = e(0) = 1 — k' and A = e{iT) = 1 + k' we obtain 

limp|fl±4^ = (l + k'r/-. (58) 

The thermodynamic limit of the double products in (53) can be calculated along the 
following chain of equalities: 

lim log ■ 



" WZ=MGi) + e{e,)){e{e[) + e{e'^)) 

f E (log(£(^*) + ^(^i)) - log(e(^D + e(^ •))) 



hm Y{me',))-me,))) 



L— >oo ' 
3 



- iog(^(%) + 1 - ^^') + E iog(^(%) + 1 + k') 



^(-<A(1 - k') + 0(1 + k')) = log(l - k'2) . 



Exponentiating we get 



Y[Z=M0i) + e{e,))m) + e{9'^)) 

Combining the latter relation with (57) and (58) we obtain (52). 

In the case of the matrix elements (54) we have to use the same arguments but instead 
of the formula (56) we have to explore the relation 

(m m-1 \ \ I 1 _ k' 

Ylog{\ + e{e,))-Y^os{X + e{9',))\ = ^ log |^ + log(A + l + k') . (60) 
1=1 i=i / 

So, in both cases for the order parameter we obtain the famous result (55). 
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6 Conclusions 



Starting from Baxter's extension of the Onsager algebra we have found factorized expres- 
sions for the spin operator matrix elements between the eigenstates of Hamiltonians of 
the finite length superintegrable ZTv-symmetric chiral Potts quantum chain up to a scalar 
factor for any pair of the Onsager sectors. 

Thus the problem of finding all the matrix elements is reduced to the calculation of 
these scalar factors. Derivation of these factors will probably require an information on 
the Bethe-states of r2-model and sZ(2)-loop algebra symmetries. Further investigation in 
this direction is important since the knowledge of the explicit formulas for the matrix 
elements will be useful for the study of the correlation functions of the superintegrable 
chiral Potts chain. These scalar factors are known for the matrix elements between the 
ground states Onsager sectors and this allows to take the thermodynamics limit using the 
standard technique and to derive the order parameters of the superintegrable chiral Potts 
chain. 

In the case N = 2 (the Ising quantum chain in a transverse field) these scalar factors 
can be found for all possible pairs of the Onsager sectors. This allows to obtain the 
factorized expression for the spin operator matrix elements in this model found in [14] in 
the framework of the method of Separation of Variables. We will address this problem in 
our forthcoming paper. 

Note added: When our paper was ready for publication we learned about Baxter's 
new paper [35], which is a continuation of his previous works [12, 13]. In this paper it is 
shown that the determinant representation (7.11) from [12] implies the factorized formulas 
for {Air ') , which in the thermodynamic limit give the order parameters of the SCPC. 

(2) 2 

Baxter's formulas for [Mi- ) coincide with our results (53) and (54). 
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Appendix A. Proof of the summation formulas 

In the case of even ^ we have the following summation formula: 



We start from the left-hand side. It does not have poles at Za = Zb- Indeed, such a pole 
may arise when a € I, b & I. But to each such summand there corresponds summand 




(_l)M/^-2)/8(^_^)/./2 



Uci^c + u){Zc + V) Uc<si^s + Zc) ■ 



(61) 
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with / replaced by / = {6} U /\{a}. A straightforward analysis shows that the sum of 
these two summands does not contain a pole at Za = Zb- Therefore the left-hand side can 
be presented as 

Pfj,{u,v;zi,Z2,. ■ ■ ,Zf,) 

where -P^(m, v; zi, Z2, ■ ■ ■ , Zp) is a polynomial of corresponding variables. Denote the degree 
of the rational function in the left hand side of (61) with respect to the variable Zc by dc 
and estimate it. From the initial sum we have dc < — /i — 1. From the factorized expression 
we have dc = deg^^ — fj, — 1. Hence deg^^ = 0. That is P^ does not depend on Zc for 
all c. 

To find Pfj,{u, v) let us analyze the pole at zi —Z2. Such pole arises if 1 € /, 2 E / 
or if 2 E /, 1 G /. Then we separate factors depending on zi and Z2 and denote by /' and 
/' the sets / and / with 1 and 2 excluded. At zi — > —Z2 the summation is reduced to the 
summation over /' reducing the calculation of the sum to the same problem for fj, — 2. It 
gives the recurrent relation 

Pp{n,v) = {-l)(>^-'y\n-v)Pp.2{u,v). 

Calculating explicitly the sum at /.i = 2 we get ^2(^*5'^) = {u — v). Therefore Pf^{u,v) = 
(_l)A'(/^-2)/8^^ — v)^/^. It proves the summation formula for even fi. 
The summation formula for odd fi 

'^I^I.(m+i)/2 ^ 

(_l)(/^+l)(/^-l)/8(^ _^ ^;)(m-1)/2 
UciZc + u){Zc + V) UcKsi^s + Zc) 

can be proved in the same way. 



Appendix B. Proof of (56) 

The aim of this Appendix is to prove (56) and (60). The proof is close to the derivation 
presented in [32]. 

We use the BAMP polynomial in e^" of degree m 

^p(.) = e--j:a.-(^^^) . (62) 
fc=o ^ ^ 

In most papers the polynomial (62) is written in terms of the variable z = ([5], similarly 
in [30] etc.). The zeros of the polynomials are known to be all simple and to come 
on the negative real axis of the variable z^ . For the following, the use of the variable u 
will be useful: so, (62) has simple zeros at u = ui,U2, ■ ■ ■ , Um, where ui,U2, ■ ■ ■ , Um can 
be chosen to have imaginary part tt/N. Each zero e^"^ of this polynomial corresponds to 
a real value of 6j, calculated through the relation 

cos0 = (l + e^")/(l-e^"). 

This relation also implies that Reu — s- — oo corresponds to 9 = 0, e{0) = 1 — k', Reu +oo 
corresponds to 9 = tt, ^(vr) = 1 + k'. 
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Figure 1: Integration contours Ci, C2, C3 used in Appendix B. 



We present the following sum as an integral along the contour Ci shown in Fig. 1: 



|^log(A + e(^,)) = £ log (3^ 



+ k' + Ay ^'p(u) 27ri 

+Ml + k' + A)/ (63) 

/ci ^p(u) 27ri 

The motivation for introducing the denominator log(l + k' + A) in the first term of r.h.s 
of (63) will be given below. The integral in the second term gives m. The other sum 
entering (56), (60), and corresponding to sector with charge P + r, is given by a similar 
relation where the second term is nlog(l + k' + A). In the case of (56) we have m = n and 
the difference of these terms disappears. In the case of (60) we have m = n + 1 and the 
difference is log(l + k' + A). 

Let us analyze the first integral from the right-hand side of (63). Its integrand is 
analytic in the domain enclosed by contours Ci and C3 except for the simple poles enclosed 
by Ci and a branch cut (where e{9) is pure imaginary) on the negative real axis along the 
segment u € [logp, 0] (enclosed by C2), where 

At n ^ 0, e{e) ^ 00. At n ^ logp, e{0) 0. 

The integral with the same integrand along the contour Ci — C3 gives 0. Indeed, the 
integral along the horizontal lines vanishes due to the periodicity of integrand under the 
transformation u ^ u + 2ti\/N . The integral along the right vertical segment vanishes 
due to 

^°<rTFTAj^°- 

This gives the motivation for introducing the denominator 1 + k' + A in the first integral 
in (63). The integral along the left vertical segment vanishes since ^ p{u) is a polynomial 
in e^" with non-zero free term and therefore '^'p{u) /'^ p{u) ~ e^" ^ as Reu — > —00. 
Therefore instead of the integration along Ci we may integrate along C3. 
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Now we shrink the contour C3 to C2, that is around the branch cut. On C2, the term 
in (62) with k = exponentially dominates as L — > 00. Indeed, the absolute value of the 
ratio of the term k = and any other term A; 7^ is the L-th power of 

|1 - cj^^e"! 

— ^ > 1 



1 - e' 

with the left-hand side of the inequality minimized at u = log p. Therefore in thermo- 
dynamic limit we can ignore the terms with /c 7^ in the polynomial (62) and replace it 

by 

e~^"[(l-e^")/(l-e")]^ 
in the integral along C2. Finally we expand the contour to C3 and get 

27ri Jc, ^\l + k' + Xj\ \ l-e-^« l-e-«yy 

The terms proportional to L dominate at L ^ 00 but they do not depend on P and so 
they disappear in (56) and (60). Therefore the left-hand side of (56) in the limit L — > 00 
becomes 

2" /c3 'vi + k'+^y 



r 



, , /l-k' + A\ r-oo+iTT/N ^ + A 

log / du = — log 

2vri ^\l + W + Xj J_^_,^/^ N ^\l + W + X 



where we used the fact that the integral along the horizontal lines vanishes due to the 
periodicity of the integrand, and along the right vertical line due to vanishing the inte- 
grand. Only the left vertical line gives a constant contribution. Recall that Ren —00 
corresponds to 6 = 0, e(0) = 1 — k'. In the case of (60) we have the same result but with 
the additional term log(l + k' + A) as was explained above. 



Appendix C. Numerical verification 

Before starting our research we implemented a numerical verification of the finite-size 
formula (53), and compared this to the analytical formula (55) for the square of the order 
parameter. We approximate the angles 9i and 0,- for L — > 00 as [34] 



sin^ {Ki + vr/iV) - sin^ {Ki 
2NI + 2Q-N 



cos 9] = TT TT , (64) 

sin 



Ki = vr , /=1,2, ...,m. 

' 2NL ' ' ' > 

We fix k' = 0.8 and obtain: 
AT = 3, L = 300, r = 1: 

Numerical calculation of (53) with (64): |q=i < S >q=2 P = 0.796894034089816 
Analytic formula: (1 - k'^)^/^ = 0.796893997784373 
Error: 3.63054 * 10"^ 

N = 5, L = 300, r = 1: 

Numerical calculation: |q=i < S >q=2 P = 0.8491969001045062 
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Analytic formula: (1 - Vl'^fl'^^ = 0.8491969000417509 
Error: 6.27554 * 10"^^ 

N = 5, L = 300, r = 2: 

Numerical calculation: |q=i < 5 >q=3 ^ = 0.7825509089016481 
Analytic formula: (1 - k'^)6/25 ^ 0.7825509088837955 
Error: 1.07617 * 10"^^ 

These results of the numerical verification encouraged us to look for an analytical proof 
of (53). 
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